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An optimization procedure for designing a 
ceramic tensile creep specimen to mini- 
mize stress concentration is carried out 
using a finite element method. The effect 
of pin loading and the specimen geometry 
are considered in the stress distribution 
calculations. A growing contact zone 
between the pin and the specimen has 
been incorporated into the problem solution 
scheme as the load is increased to its full 
value. The optimization procedures are 
performed for the specimen, and all 
design variables including pinhole location 
and pinhole diameter, head width, neck 
radius, and gauge length are determined 
based on a set of constraints imposed on 
the problem. In addition, for the purpose of 



assessing the possibility of delayed fail- 
ure outside the gage section, power-law 
creep in the tensile specimen is consid- 
ered in the analysis. Using a particular 
grade of advanced ceramics as an 
example, it is found that if the specimen is 
not designed properly, significant creep 
deformation and stress redistribution may 
occur in the head of the specimen result- 
ing in undesirable (delayed) head failure of 
the specimen during the creep test. 

Key words: ceramics; finite element 
method; optimum design; power-law 
creep; specimen design; stress analysis; 
stress relaxation; 

Accepted: February 15, 1996 



1. Introduction 



Because of their high temperature strength, oxidation 
and corrosion resistance, and superior wear perfor- 
mance, advanced ceramics (e.g., sihcon nitride, sihcon 
carbide, alumina, etc.) are promising materials for use 
in high temperature, load-bearing applications such as 
turbine engines, heat exchangers and waste incinerators. 
Since these devices are designed for an extended service 
life, their creep resistance and long-term reliability 
under service must be ascertained. Accordingly, stress 
and temperature dependent creep data are needed for 
designing a structural component for long-term usage at 
elevated temperatures. 

In the ceramics community, flexural creep testing 
using three-point or four-point bending provides an easy 
way of collecting creep data, as this type of loading does 
not create problems in alignment and gripping. How- 



ever, to obtain reliable creep data from these bending 
specimens which contain non-uniform stress fields, a 
proper data interpretation method must be implemented. 
Assuming the material follows a power-law creep 
characteristic, a deconvolution method is available for 
sorting out the asymmetric creep behavior from the 
bending data to yield both tensile and compressive creep 
data [1-7]. 

In order to bypass the cumbersome task of decon- 
volution and data interpretation, a direct, simple tensile 
or compressive creep test may be adopted, even though 
critical issues such as alignment, fixturing, and gripping 
must be faced and dealt with. In the case of uniaxial 
tensile creep, there are a variety of testing methods that 
differ in terms of the design of the specimen geometry, 
fixturing, gripping, and the displacement measurement 
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method (see, for example, Refs.[8-14]). All techniques 
have their advantages and shortcomings. We take the 
experimental set-up at NIST [7, 10, 15] as an example 
(see Fig. 1). The adopted geometry of a typical tensile 
creep specimen is of the coupon type consisting of two 
heads, a gauge section, and two ramps (or necks) 
connecting the two entities. A central circular hole is 
drilled precisely at the centerline in both heads from 
where the loading pins can run through to assure a good 
alignment. A constant load is applied at the loading pins 
and the separation of the two flanges attached at the 
gauge section is continually monitored using a laser 
beam for creep strain measurements [10]. In an ideal 
case the specimen should rupture somewhere within the 
gauge section at the end of the experiment or remain 
intact in an interrupt creep test as the nominal stress is 
highest, and is uniform throughout the whole gauge 
section. However, in many cases, the specimen breaks 
prematurely during the initial loading stage. Machining 
damage is one possible cause of the undesirable failure 
during fatigue testing [16], but more often premature 
failure is observed from the pinhole in the head or in the 
neck, presumably due to elastic or creep stress concen- 
tration. Since a hole or a neck is a well-known stress 
concentrator, a proper design of the specimen to reduce 
stress concentration to a minimum becomes an impor- 
tant task, so that the probability of rupture is enhanced 
in the gauge section in relation to other places (i.e., neck 
or head sections). 



The objective of the present study was to launch an 
optimum design of tensile creep specimens of the type 
sketched in Fig. 1, so that the chances of undesirable 
failure can be reduced to a minimum. Owing to the 
complexity of the specimen geometry, finite-element 
analysis offers a viable vehicle to undertake the opti- 
mum design task. In this paper, we use a commercial 
finite element package, ANSYS , ^ to perform the 
optimization analysis for the design of this tensile 
specimen. In order to obtain the stress distribution 
around the loading pinhole, contact elements were incor- 
porated into the model for stress analysis. The creep 
analysis was carried out to investigate the stress redistri- 
bution and the creep deformation in high strain regions. 

2. Original Designs 

Before we initiate the task of optimization, we briefly 
review the original design based on many years of expe- 
rience accumulated at NIST. Figure 1 describes the 
original design of the tensile creep specimen geometry 
[7, 10, 15]. There are three different designs of speci- 
men geometry in terms of total length TL ; pin-hole size 
D\ neck radius R\ head width W\ and gage length GL^ 
etc. as tabulated in Table 1. For example. Specimen 1 
has the following dimensions: 30.00 mm total length, 
2.00 mm thickness, 10.00 mm gauge length, together 
with a 7.00 mm head width, 2.44 mm pinhole diameter, 
and 2.50 mm neck radius. A constant load is applied to 



POINT 1 POINT 2 POINT 3 POINT 4 



POINT 5 




Fig. 1. Initial design of the tensile creep specimen. Point 1: Head stress concentrator; 
Point 2: Hole stress concentrator; Point 3: Neck stress concentrator; Point 4: Gage section; 
Point 5: Load point. 



Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 
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Table 1. Specimen dimensions in the initial design 



Specimen 


TL 


GL 


D 


R LI 


W 


GS = BXT 


No. 


(mm) 


(mm) 


(mm) 


(mm) (mm) 


(mm) 


(mm) 


1 


30 


10 


2.44 


2.5 4.0 


7.0 


2.0X2.0 


2 


50 


20 


4.76 


2.5 6.25 


12.5 


2.5X2.0 


3 


76 


19 


4.78 


19.0 8.25 


15.9 


2.5X2.5 



Notations: TL = total length; GL = gage length; D = hole diameter; 
R = neck radius; LI = hole position; W = head width; GS = gage sec- 
tion (width X thickness, BXT); see Fig 1. 



the specimen from the pinhole (pin-loaded) via a 
shackle or clevis with a 3.05 mm opening and a pin rod 
of 2.44 mm diameter slide fit with a tolerance of 
±0.025 mm. The whole assembly including the speci- 
men and loading fixture is enclosed in a furnace for high 
temperature testing. Windows in the furnace allow the 
laser beam to go through for the purpose of monitoring 
the creep displacement. Because of the experimental 
set-up, the working space is somewhat limited; accord- 
ingly the total length of the specimen is set at a fixed 
value of 30 mm, 50 mm, or 76 mm. The current design 
has been used to produce tensile creep data for silicon 
carbide [12], silicon nitride [17], and other ceramic 
materials [18]. Sometimes the specimen breaks at the 
neck or head region at the initial loading stage. Figure 2 
shows one example of such a failed sample. In order to 
prevent premature failure of the specimen, redesigning 
the specimen is clearly warranted. 








llflMillM 



3 29 30 31 32 33 34 ^ 




^A --t:^?*^ 



(1989-90 vintage) Inc. for the specimen and sintered 
silicon carbide with a trade name of Hexolloy produced 
by Norton Ceramics Co. for the loading pin rod. Their 
mechanical properties at ambient temperature, as sup- 
plied by the vendors, are as follows: 

For Si3N4: modulus of elasticity, E = 350 GPa 
Poisson's ratio, v = 0.24 
ultimate tensile strength, ^u = 900 MPa 
with a Weibull modulus =15. 

For SiC: modulus of elasticity, E = 400 GPa 
Poisson's ratio, v = 0.24 
ultimate tensile strength, ^u = 500 MPa 
with a Weibull modulus = 17. 

These materials were chosen from a round robin testing 
program in which NIST participated. 

In the present paper, we will present a two-part 
design analysis. In the first part, we analyze the stress 
and strain fields around the pin hole caused by the 
contact between the pin and the specimen. This will help 
us understand the cause of premature failure for many 
specimens of original design, and plan our design path 
in the second part of the analysis. In the second part, we 
systematically alter the geometry to minimize the stress/ 
strain at the critical spots, thereby yielding the final 
design for the optimized geometry. 



3. Optimum Design Procedure 
3.1 Initiation of Optimum Design 

As stated, throughout the entire design procedure, the 
total length of the specimen will be fixed at either 
30 mm, 50 mm or 76 mm depending on the specimen 
class. Also, the opening of the clevis will be fixed at 
3.05 mm. These two conditions are regarded as given, 
and are the only two constraints imposed on the prob- 
lem. All other parameters will be varied to find the 
optimum dimensions. In the optimization scheme, we 
will take the current design parameters shown in Fig. 1 , 
except the pinhole diameter, as the initial values in the 
iterative optimization loop. Throughout the entire creep 
test, since the applied load is kept constant, the nominal 
tensile stress in the gauge section will always be main- 
tained at 100 MPa. 



Fig. 2. Photograph of a broken alumina tensile creep specimen after 
100 h test duration. 



From the current NIST tensile creep test program, the 
materials selected for study are a grade of commercially 
available HIPed silicon nitride with a trade name of 
PY-6 produced by GTE Laboratories, Waltham, MA 



3.2 Design Criteria 

In the optimum calculations, we set the following 
design criteria to be met in each iterative loop: 

(a) Smallest stress concentration in the neck area (i.e., 
the transition region from the neck of the flange to 
the gauge section). 
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(b) Tensile stresses as high as possible in the loading 
pin rod, but without failures. 

(c) Smallest critical stresses in the head. 

(d) Longest possible gauge length. 

4. Design Analysis 

4.1 Pin-Specimen Contact Analysis 

When we investigate the stress distribution on the 
specimen, particularly around the pinhole in the head, 
the effect of the contact surface between the pin and the 
pinhole must be considered because, for a given applied 
load, the traction forces will depend on contact area and 
the local stresses will be influenced by the applied trac- 
tion at the pinhole. In the ANSYS software, a 2-D Point 
to Surface Contact Element named CONTAC48 [19] 
intended for general contact analysis is available. The 
contact element contains three nodes, one at the contact 
surface and the other two at the target surface, so that 
either a 4-node or 3-node structural element can be used 
to model the contacted region. 

Owing to the symmetry conditions, only one quarter 
of the specimen geometry needs be considered in the 
finite element modeling. The overall 2-D finite element 
mesh for the quarter-specimen is shown in Fig. 3, taking 
the 50 mm originally designed specimen as an example. 
For loading conditions, we consider a point load applied 
at the right side of the pin as indicated in Fig. 3. 

4.2 Creep Analysis 

Since the tensile specimens are intended for use at 
elevated temperatures, a creep stress and strain analysis 
is needed to evaluate the specimen behavior during 
loading. The need for optimization analysis can be seen 
from Fig. 2 where a creep specimen failed prematurely 
in less than 100 h of the creep test. An undesirable 
failure has occurred from the pinhole, which is mainly 
induced by a significant creep deformation in the speci- 
men head. Therefore we must pay more attention to the 
creep behavior in this region. 



The finite element package ANSYS has a feature for 
creep calculations for most engineering materials [21] 
and various formulations describing creep law could be 
directly used in the creep analysis. Here we choose an 
Arrhenius type of equation in which the values of the 
coefficients are provided for hot-pressed silicon nitride 
from the literature [17]. 



^s=Ao(cr/cro)"e RT 



(1) 



where g^ is the steady state creep strain rate, ao is the 
reference stress, Ao and n are materials constants, a is 
the applied stress, Q is the apparent activation energy for 
creep, R is the universal gas constant, and T is the abso- 
lute temperature. For the purpose of computation, we 
consider the isothermal conditions at r= 1600 K in 
which Q= 1310 kJ/mol for the PY-6 material according 
to Krause and Wiederhorn [17] . Then, Eq.(l) reduces to 
the following simple power-law form: 



.,=A(a/aor 



(2) 



This is the well-known Norton's equation where the 
material constants A = 5.0X10"^^ h"^ and the stress 
exponent n = 8.4, if g^ is expressed in the units of (1/h) 
and (To is 1 MPa. 

In order to focus our emphasis on the pin hole area, 
we refine the local mesh in the maximum stress region. 
We choose a special type of element available in ANSYS 
with analyzing creep capabilities to handle the time- 
dependent creep strain and stress redistribution compu- 
tations. 

4.3 Steps in an Iteration 

In the ANSYS software, a design optimization routine 
is available. The routine is capable of performing a series 
of analysis-evaluation-modification cycles until all 
specified design criteria are met [20] . The procedure of 
design optimization consists of the following six main 
steps: 




Fig. 3. Two-dimensional finite element mesh for the quarter-specimen. 
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(a) Initialize the design variables In this paper, we 
performed the optimization for the tensile specimen 
geometry. As was shown in Fig. 1, there are four design 
variables: the gage length GL, the radius of curvature of 
the neck R , the width of the head W, and the location of 
the hole center L 1 . The initial values of these design 
variables, which represent the starting point in the 
design loop, will be later modified by the ANSYS 
optimizer. 

(b) Build the model parametrically In this step, we 
construct the model in terms of the design variable 
parameters. An 8-node structural plane element named 
PLANE82 [22] is used for the specimen design. 

(c) Acquire the solution In this step, we first define 
the analysis type as the linear static analysis and set the 
appropriate applied load level necessary to make the 
uniform tensile stress of 1 00 MPa in the gage section of 
the specimen. Then we seek the finite element solutions 
for each intermediate design. 

(d) Retrieve the results parametrically and set the 
state variables and objective function parameters In 
the case of the specimen design, there are three state 
variables: namely the hole stress, neck stress, and head 
stress . 

(e) Declare optimization variables and begin the 
optimization process The optimization paths for the 
specimen design are described in the next section. 

(f) Review and verify the results We review and 
verify the results of the optimization run by plotting the 
graph of state variable versus design parameter, so that 
we can track how a variable changed from loop to loop. 

4.4 Optimization Paths 

There are at least four independent design variables 
that must be determined at the end of the optimization 
procedure for the specimen design: head width W, gage 
length GL, head length LI, and neck radius R. The 
optimization procedure is divided into three parts with 
three different objective functions and dominate 
parameters. In the first part, designing by pinhole stress, 
we use Wl = W-D as an objective function dominated 
by maximum pinhole stress at the edge of the hole. In 
the second part, designing by head stress, we use 
Lii = Ll-D/2 as an objective function dominated by 
head stress. In the third part, designing by neck stress, 
we use L2 = TL-{GL^-R+L\) as an objective function 
dominated by neck stress. 

5. Results 

5.1 Elastic Stresses in the Specimen 

Analysis results of the stress distributions for the 
quarter-specimen in the case of the original design for a 



50 mm specimen are shown in Figs. 4a and 4b. Figure 
4a is the principal stress contour plot based on contact 
analysis. The solution of the equivalent stress distribu- 
tion is given in Figure 4b, where the equivalent stress is 
defined hy a^=v (3o-'ijCr\j/2X where a\j is the devia- 
toric stress tensor, and the repeated index denotes sum- 
mation. As can be seen in the figures, a maximum 
principal stress is found on the edge of the hole 
which is a stress concentrator (see Fig. 4a). There are 
two additional stress concentration zones; one is in the 
neck area located at the intersection of the straight 
gage section and the curved neck region, and the other 
is located in the head at the end of the specimen 
along the centerline. Moreover, there exists a maximum 
equivalent stress at the pin and pinhole edge in 
Fig. 4b. 

5.2 Pin and Pinhole Dimensions 

Since it is well-known that increasing the pin diame- 
ter will result in a monotonically decreasing stress on 
the pin rod, it follows that the sizes of both the pin rod 
and the slightly larger pinhole should be as large as 
possible to the extent imposed by the finite physical 
dimension. This result is applicable to the general case 
of ceramic materials. For the case of a silicon carbide 
rod and silicon nitride specimen, the general rule of 
thumb is that the normalized pin diameter is approxi- 
mately 0.1 (i.e., D/TL = OA) of the total specimen 
length. For example, in the case of the pin rod designed 
originally, D = 4.76 mm for the 50 mm specimens. 
Hereafter we will use this constraint to perform the 
optimization task for the 50 mm specimens, as an illus- 
trative example. 

5.3 Specimen Geometry 

(a) Design by hole stress Taking a series of hole stresses 
against W, we obtain Fig. 5a. If an allowable value of 
hole stress is specified for the specimen design, we can 
then determine the W value. In the same way, we get 
another group of values of hole stresses from Fig. 5b for 
a range of hole locations and determine the optimum 
value of hole position L 1 . By the same token. Fig. 5c is 
a plot of hole stress for a range of gage lengths, and the 
best gage length can then be determined from this plot. 
It turns out that the stresses are monotonically increas- 
ing for increasing gage lengths. For the sake of minimiz- 
ing stress, the smallest gage length should be used. 
However, to fully utilize the material, the longest gage 
length should be adopted. Clearly a balance must be 
struck between these two opposing ''forces." Similarly, 
Fig. 5d plots the hole stresses with increasing neck 
radius. 
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Fig. 4a. Stress distribution from the contact analysis: the first principal stress. 
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Fig. 4b. Stress distribution from the contact analysis: the equivalent stress. 
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Fig. 5a. Hole stress of the specimen as a function of head width W 



Fig. 5b. Hole stress of the specimen as a function of hole location LI. 
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Fig. 5c. Hole stress of the specimen as a function of gage lengtli GL 



Fig. 6a. Head stress as a function of liole position L 1. 
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Fig. 5d. Hole stress of the specimen as a function of neck radius R. 



Fig. 6b. Head stress as a function of head width W. 



(b) Design by head stress For a given value of pinhole 
diameter D = 4.76 mm. Fig. 6a shows the relationship 
between head stress and hole position LI. The head 
length 2XL1 can be easily chosen from Fig. 6a, if an 
allowable value of the head stress is imposed for the 
design. Similarly, Fig. 6b plots the head stress vs head 
width W for the given D , and the optimum value of W 
can thus be chosen. 

(c) Design by neck stress From Fig. 7, we find that the 
neck stress is strongly dependent upon the value of the 
neck radius R , so that if the allowable value of the neck 
stress is given in the design, the neck radius can be 
determined accordingly. 
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Fig. 7. Neck stress vs neck radius R. 
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5.4 Allowable Stresses 

A series of criteria must be given for the specimen 
design, including the allowable value of hole stress a\oh 
on the edge of the pinhole, the allowable value of head 
stress cr'^head, and also the allowable value of neck stress 
c'^neck- The design criteria are expressed as follows: 



f^head — ^ head 
f^neck — ^ neck 
f^neck — ^ neck 



(3) 
(4) 
(5) 



where a"" is the allowable stress. In the optimum design, 
we set the following criteria after taking the allowable 
stress values into account for the selected silicon nitride: 



0-\ole = O-'neck = ^-'head = 105 MPa. 



(6) 



We then obtain the following set of design parameters 
for a 50 mm total length specimen after a series of 
iterations: 



D =4.76 mm 
W = 12.5 mm 
R =9.0 mm 
LI =7.1 mm 
GL = 12.0 mm 



(7) 



The same procedures have been carried out on the other 
specimens, i.e., a 30 mm total length specimen and a 
76 mm total length specimen. The final results are given 
in Table 2 in nondimensional ratios which are formed 
via normalization to the total specimen length TL. 

Table 2. Sizes for final optimum design 



D 


GL 


R 


W 


LI 


B 


0.10 


0.25 


0.20 


0.25 


0.14 


0.05 



^ Dimension normalized to total length TL of a specimen. 



5,5 Creep Behavior 

The stress redistribution contour for the quarter-spec- 
imen due to creep is displayed in Fig. 8. At the early 
stage of creep, the stress near the hole edge is high but 
drops dramatically with time. The stress relaxation 



behavior due to creep is shown in Fig. 9. At /^ = 5 h, the 
initially designed specimen (Fig. 9a) for example, the 
stress at a distance 0.27 mm from the hole is observed to 
have the highest stress (106 MPa), then drops to 84 MPa 
at t = 200 h. For a finally designed specimen, however, 
the stress level at the same location (0.27 mm) has 
diminished to about 80 MPa, and the highest stress 
(104 MPa) at the edge of the hole reduces to 76 MPa at 
t = 200 h (Fig. 9b). Overall, we see a gradual relaxation 
of the stress concentration around the hole for both 
designs of specimens with reduced stresses for the final 
design (see Figs. 10a and 10b). 

The stress redistribution in the critical area (i.e., at the 
hole stress concentration zone, gage length nominal 
stress section, neck stress concentration zone, head 
stress concentration zone, and the load-point applied 
stress area) due to creep are displayed in Figs. 1 la and 
lib, respectively, for the initial and final designs of the 
specimen. As can be seen for both cases, relaxation 
takes place at the hole edge and neck root but increases 
slightly at the head and load-point, whereas it remains at 
an approximately constant level within the gage section. 
Again, we see that the creep stresses in the final design 
have been dramatically reduced. As opposed to the 
elastic case, the stress peak does not occur on the edge 
of the hole, after the stress redistribution is modified by 
creep. From Fig. 8, we also observe that a bigger 
equivalent stress (see Fig. 8b with t= h. Fig. 8d with 
r = 50 h, and Fig. 8f with t = 200 h) is exactly located at 
the place where creep fracture may occur in the 
specimen. The analysis seems to suggest that creep 
failure follows a maximum strain criterion. Hence, 
if a design against early failure due to creep is to 
be implemented, then the stress redistribution and 
the total strain accumulation must be evaluated (see 
Figs. 10 to 12). 

The creep strain solutions in the critical areas are 
given in Fig. 12 for both cases. As can be seen, after 
about 100 h the creep rates (i.e., the slopes) around the 
hole are almost identical. However, for the initial design, 
the accumulated strains at the hole edge are very high 
within the first 80 h (Fig. 12a). If the failure criterion of 
this material is such that failure will occur when the 
local strain exceeds the critical strain (or creep ductil- 
ity), then the cracking from the pinhole will be observed 
at the location corresponding to this element. However, 
from Fig. 12b we see that the total creep strain at the 
hole edge has been reduced from 0.044 % to 0.015 %. 
Thus, based on the critical strain failure criterion, the 
optimized specimen should not suffer premature failure 
during creep test. 
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Fig. 8a. Creep stress redistribution in contour plot: principal stress at time = 0. 
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Fig. 8b. Creep stress redistribution in contour plot: equivalent stress at time = 0. 
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Fig. 8c. Creep stress redistribution in contour plot: principal stress at time = 50 h. 
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Fig. 8d. Creep stress redistribution in contour plot: equivalent stress at time = 50 h. 
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Fig. 8e. Creep stress redistribution in contour plot: principal stress at time = 200 h. 
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Fig. 8f. Creep stress redistribution in contour plot: equivalent stress at time = 200 h. 



1B0 



□ 


D - mm 


o 


D ■0.27 mm 


o 


D - . ^2 mm 


* 


□ I. 0.57 mm 


K 


D = . 72 mm 


M 


Dw0.53mrtl 




50 100 150 

TIME (HOURS) 



200 



Fig. 9a. Stress relaxation at the region near pinhole as a function of 
time for initial specimen design. 
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Fig. 9b. Stress relaxation at the region near pinhole as a function of 
time for final specimen design. 
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Fig. 11a. Creep stresses experienced at various critical points as a 
function of time for initial specimen design. 
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Fig. 10b. Time-dependent stress distributions in the near-hole zone 
for final specimen design. 
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Fig. lib. Creep stresses experienced at various critical points as a 
function of time for final specimen design. 
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Fig. 12a. Creep strains experienced at the corresponding locations 
shown in Fig. 1 for initial specimen design. 
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Fig. 12b. Creep strains experienced at the corresponding locations 
shown in Fig. 1 for final specimen design. 



6. Discussion 

6.1 Optimum Design Procedure 

In the process of optimization, we have imposed the 
constraint that the total length of the specimen be fixed 
at either 30 mm, 50 mm, or 76 mm (see Table 1). The 
length of the specimen usually is dictated by the size of 
the billet, fixture design, and other environmental 
factors. Since the current optimization scheme is based 
on elasticity analysis, the geometry of the specimen 
under investigation is scalable, and nondimensional 
geometric parameters normalized by the total length 
may be used. For instance, instead of searching for the 
optimum head width V^ , the nondimensional width 
w = W/TL can be used. Table 2 lists the results of 
optimization on a group of geometric parameters. These 
results can be scaled-up to acquire the optimum dimen- 
sions not only for the 30 mm, 50 mm, and 76 mm 
specimens, but also for other specimens with different 
lengths. 

6.2 Double-Reduction Design 

Sometimes a situation may arise such that after the 
optimum design is implemented the total load applied 
through the pin is so high that the pin fractures first. One 
possible remedy is to use higher strength pins, or alter- 
natively, the total load could be reduced. In the latter 
case, to reduce the total load while maintaining the 
same level of tensile stress in the gage section, a double- 
reduction design (to reduce the gage cross-section in the 
z -direction) may have to be implemented. This means 
that a three-dimensional finite element model must be 
constructed in order to perform the optimum design. A 
3-D analysis is possible, but is more labor-intensive and 
time-consuming. We regard this as outside the scope of 
the present paper and thus we will not pursue this 
subject further. 

6.3 Creep Asymmetry 

It is known that ceramic materials usually creep less 
when subjected to compressive stress than tensile stress. 
This asymmetric creep behavior certainly applies to 
many grades of silicon nitrides. There are a few excep- 
tions such as hot pressed alumina which shows symmet- 
ric behavior [23]. In the present paper, the creep proper- 
ties had been assumed to be symmetrical for the ease of 
computations. Our results showed that for the present 
material being analyzed, creep symmetry seems to give 
rise to delayed head failure. This conclusion is consistent 
with the observation of the resulting delayed failure of 
alumina presented in Fig. 2. Moreover, no delay failure 
has been observed in the case of silicon nitride and 
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silicon carbide so far which had been known to exhibit 
asymmetric creep behavior. Thus, it is likely that creep 
asymmetry may have played a beneficial role in prevent- 
ing delayed head failure of the creep specimens. Proof 
of this contention should be an interesting task for future 
research. 

6.4 Failure Mechanisms 

Upon completion of the optimum design, specimens 
were prepared based on the dimensions obtained except 
the 76 mm specimen [24]^, and the tensile creep testing 
on these specimens was carried out. The results showed 
that the failure rate outside the gage section has been 
reduced significantly. In the rare cases where the 
undesirable failure did occur shortly after applying a 
load, it is always an elastic failure due to inherently 
existing natural flaws in the material as received, which 
is normal for ceramic materials. In cases where delayed 
failure occurs, the fracture surface showed different 
characteristics, indicating that creep damage starts to 
accumulate. The failure in this case seems to follow the 
critical strain criterion, as creep damage builds-up to the 
critical strain level. 

7. Summary 

We have performed an optimum design of a tensile 
creep specimen in a plate-type geometry using the 
commercially available finite element code ANSYS as 
the design tool. We set the initial conditions based on the 
current design being used in the NIST program. Two 
constraints were imposed on the optimum design calcu- 
lations: (1) the total length of the specimen is set at a 
fixed length (say, 30 mm, 50 mm, and 76 mm); and 
(2) the opening of the clevis is set at 3.05 mm. Adopting 
an allowable stress of 105 MPa, the results of the 
iterative computations yield the following values regard- 
ing the final design of the specimen: W/TX = 0.25, 
DITL = 0.10, L l/TL = 0.14, R/TL = 0.20, and 
GLITL = 0.25. It is expected that this optimum design 
should give the smallest probability of premature failure 
as frequently encountered in the current design. 

In addition, we also performed a complete time- 
dependent creep analysis regarding the evolution of 
stress and strain as functions of time-history. The results 
indicated that the peak stresses and stains occurred 
at locations different from those predicted by elastic 
analysis. 



Since the completion of the manuscript, the dimensions of the 
30 mm specimen have been altered shghtly due to enlargement of the 
head size: W^ll.Q (7.5) mm; D = 3.57 (3.00) mm; /? = 5 (6) mm; 
LI = 4 (4.2) mm; GL = 8(7.5)mm, where the numbers within the 
parentheses designate optimum dimensions. 



The predicted critical location coincides with the 
location where rupture actually occurred in an alumina 
specimen that failed prematurely before 100 h. The 
creep analysis seems to suggest that creep failure follows 
a maximum strain criterion. Comparisons of creep 
stresses between as-received and optimized specimens 
showed reduced severity, thereby lowering the probabil- 
ity of premature failure during creep test for the final 
designed specimens. 
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